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Abstract 

In our five-dimensional cosmological model, we investigate the role of a Lorentz vio- 
lating vector "aether" field on the moduli stabilization mechanism. We consider the 
case of a space-like ^ther field on a compact circle with Maxwell-type kinetic term. 
The Casimir energy of certain combinations of massless and massive bulk fields gen- 
erates a stabilizing potential for the radius of the compact direction while driving the 
accelerated expansion in the non-compact directions. It is shown that the «ther field 
can reduce the influence of the Casimir force and slow down the oscillation of the 
radion field. This property proves crucial to the stability of the extra dimension in 
the universe where non-relativistic matter is present. We speculate that this scenario 
might reveal a hidden connection between the dimensionality of spacetime and the 
spontaneous breaking of Lorentz symmetry. 
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1 Introduction 



Although there have been a lot of progresses on constructing phenomenologically viable 
models based on theories with extra spatial dimensions, some fundamental questions have 
not been completely solved. One of them is the moduli stabilization problem. The size and 
shape of compact space described by dynamical moduli fields have to be fixed in order to 
avoid any conflict with astronomical observations. In addition to these problems, we also face 
the challenges from cosmology in explaining the accelerated expansion of the universe. One 
possible solution for these problems may involve arguments based on anthropic principle. 
However, the search for an alternative solution is still going on, for example in [1]. 

Recently, it was suggested that Casimir energy from various field fluctuations in compact 
extra dimensions could play a crucial role in addressing these significant problems [21 E]- 
Greene and Levin [3] argued that if the total Casimir energy is properly chosen, then it is 
possible, at least in the case of vacuum dominated universe, to stabilize the size of the extra 
dimensions and drive the accelerated expansion of the three non-compacted directions in 
which the Casimir energy plays the role of dark energy. The authors in [1] employed the 
calculation of Casimir energy in the non-trivial space, M^+^ x and demonstrated that 
the shape of extra dimensions can also be stabilized by the same mechanism. Interestingly, 
predictions in this scenario such as radius of extra dimensions and quantum gravity scale 
in the bulk are in agreement with those from the large extra dimensions or ADD scenario 
[21 E]- Hence, the moduli stabilization problem, dark energy problem and the hierarchy 
problem may possibly be explained in the single unified framework. However, as it was 
pointed out in [3], there are some crucial obstructions to realizing a phenomenological viable 
version of this scenario. One of them is that the extra dimension fails to stabilize if we 
include contribution from matter contents. During the matter dominant epoch, energy 
density of non-relativistic matter was the dominant contribution in the effective potential 
of the moduli fields and washed away the minimum of the effective potential. Although 
the minimum reappears in the vacuum dominated epoch, the moduli (i.e. radion field) has 
already passed the dynamical stable fixed point. This causes the extra dimension to expand 
and contradicts with our observation. Thus, it would be interesting to investigate whether 
this technical problem could be solved. 

In this paper, we propose the new stabilization mechanism based on the Casimir energy 
and the existence of the Lorentz violating "^ther" field in the compact direction. Starting 
with the simplest model with one extra dimension where the space-like aether field lives in 
the compact circle similar to the model considered in [H [8], we claim that non- vanishing 
vacuum expectation value (vev) of the aether field would affect the dynamical equation of 
background moduli field. It can reduce the gradient of the radion's potential and slows 
down the oscillation frequency. This ensures stability of extra dimension although there is 
non-relativistic (dark) matter in the universe. As in the previous works, the Casimir energy 
of massless and massive fields embedded in five- dimensional space play a role of dark energy 
and drive the expansion of non-compact space as expected. Note that the effect of a time-like 
ffither field on slowing down the expansion rate of the universe was pointed out in [9j. The 
effects of the aether fields on cosmological observable was studied in [lOj. The authors in 
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[TT| \n\ also studied the role of aether field on the stability of the extra dimension in the 
context of braneword scenario but in a different aspect. 

It is important to state that we are aware of the stability issue for space-like ^ther field 
p!3| m] which may cause difficulty in the construction of a more realistic model of this 
scenario. However the interplay between the ^ther field and the dynamical moduli field in 
our model may shed some light on the connection between the dimensionality of spacetime 
and the violation of Lorentz symmetry. Perhaps nature allows us to observe only the large 
three-dimensional space that preserves Lorentz symmetry but conceals the Lorentz violating 
directions in the compact space. 

This paper is organized as follows. In section [21 we start by reviewing the aether model 
in 5-dimensional spacetime. In section [31 we derive cosmological equations of motion in 5- 
dimensional spacetime with aether field and write down effective 4-dimensional equations of 
motion in the radion picture. Then we review the calculation of Casimir energy and extend 
to the case involving interaction between aether field and bulk fields in section [H In section 
[5l we investigate the role of aether field in the context of stabilization of the extra dimension 
both in the vacuum dominated universe and in the universe with non-relativistic matter. 
Finally we summarize our results in section [21 

2 ^Ether Field and Its Interactions 

We start by considering a 5-dimensional fiat spacetime with coordinates = {x^,y) where 
fi = 0, . . . , 3 and with mostly plus metric signature. We assume that the fifth direction 
is compactified on a circle. Now we consider a toy model in which Lorentz symmetry is 
spontaneously broken by the ^ther field u"" i.e. a vector field with a non-vanishing expec- 
tation value. Most of the aether models contain kinetic term that makes their Hamiltonian 
unbounded from below and their stability is a subtle issue [TB] . Here we consider the action 
with Maxwell-type kinetic term [8] 

S = j (fx^g ( - -^VabV''' - Kuau'' - + J] . (1) 

Here Vab = ^aUb — ^bUa has a familiar form to the field strength tensor of electromagnetism. 
However, the ^ther field is not related to the electromagnetic vector field A°- and its 
dynamics does not respect U{1) gauge symmetry. In contrast, the second term in the above 
action enforces the aether field to have a constant norm 

u'^Ua = V", (2) 

where A acts as a Lagrange multiplier and we take f ^ > 0. In our unit m° has dimension of 
mass^/^. The sum Ci in ([1]) represent various interaction terms which couple the aether field 
to matter fields that we will discuss later in this section. If we neglect the interaction terms 
for the moment, the equations of motion for the aether field u'^ can be written as 

VaV^ + v~'^u\^V dV^'^ = 0. (3) 
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Any solutions for which Vab = will solve the equation of motion In order to preserve 
Lorentz invariance in the 4-dimensional non-compact space, we choose the background so- 
lution such that the ^ther is a space-like vector field which has non-vanishing components 
along the extra fifth dimension, 

«'^ = (0,0,0,0,1;). (4) 

It is important to note that there is a subtle stability issue here. Although our aim is to 
investigate the role of the aether field on the stability of the extra dimension, the model of 
space-like aether field with Maxwell-type kinetic term itself is unstable [13] . However, for our 
purpose, we can consider it as a toy model and assume that there is some mechanism which 
stabihzes the aether field. 

The energy-momentum tensor of the aether field Tab\u takes the following form 

Tatlu = VacV,' - ^V.dV'gab + V-^UaUmc^ . (5) 

Note that properties of the aether field depend crucially on spacetime geometry. The flat 
space background solution in equation (jl]) gives Tah\u = 0. However, in curved spacetime, 
the sether field can give rise to non-vanishing energy momentum tensor for example a time- 
like aether field can produce energy density [9] while a space-like aether gives the stress 
components [14]. The case of an aether field oriented along the compact extra dimension 
was investigated in [8]. It was shown that such aether configuration can also give rise to 
non- vanishing energy momentum tensor. However, Tab\u vanishes when the extra dimension 
is stabilized. We will review this result in the next section. 

We now consider the effect of the interaction term ^ . £j in ([T]) which in general can 
include the terms corresponding to the aether field coupled to scalars, fermions and gravity. 
However, we will consider stabilization mechanism of the extra dimension involving Casimir 
energy of gravitons, bulk scalars and bulk fermions. We will ignore the bulk vector terms. 
Let us begin with the effect of the interaction of the ^ther with a real massive scalar field 
0. The Lagrangian for the scalar field with the minimal coupling term is 

^<t> = -l{da(f)Y - \m^'' - l^u^u'>da(t>db(l>^ (6) 

where /x^ is the coupling parameters with dimension of mass^/^. The corresponding equation 
of motion for the scalar field takes the form [8] 

- = -fi^^da{u''u'db(f)). (7) 

Expanding the scalar field in Fourier modes oc e*'^"^", we obtain the modified dispersion 
relation, 

-fc% = m2 + (l + aj)fc2, (8) 
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where the dimensionless parameter = v / n^f, the ratio of the aether vev to the couphng 
Next we consider the fermion terms. The Lagrangian for fermionic field with the minimal 
coupling term can be written as [8] 



where /i,^ is the fermionic coupling constant with the unit of mass^''^. In the same spirit 
as in the scalar field case, the corresponding modification of the dispersion relation for the 
fermionic case can be written as 

-kf^k^ = m^ + {l + aiykl (10) 

where the dimensionless parameter = v/fi^. The form of this equation is different from 
the analogous equation in the bosonic case: i.e. the second term on the right-handed side 
increases by instead of a^. Finally, we consider the aether field which couples non- 
minimally to gravity. This can be described by the action [S] 



d'x^/^i^j^R + agU'^u'Rab) , (11) 



where ag is the dimensionless graviton coupling constant and M^, is the Planck mass in 5 
dimensional space-time. By varying this action with respect to the metric tensor, we obtain 
the equation of motion Gab = ^''^GTab\(GC) with 

TabliGC) = ag(^R,dU'u'^gab + VcVaiUbU') + VbVciUaU') - V diu'u"^) Qab ' V cV%UaUb)^ , (12) 

where G is the 5-dimensional gravitational constant. Let us consider small fluctuation of the 
metric 

gab = Vab + hab- (13) 

Following the explanation in [8], the metric perturbation can be decomposed into 

V = V + ^V^.u, K5 = ^, (14) 

where rj^'^h^v = 0, hf^^, presents the propagating modes of the gravitational wave, $ denotes 
the Newtonian gravitational field and ^ is a component associated with the radion field 
describing the modes of the extra dimension. By setting ^ = = ^, and considering 
transverse waves, d^h\^ = 0, the gravitational equation of motion becomes 

- = ^rc^dfh,.. (15) 

2 

Let us define ctg = IQn^j^. The above equation gives the modified dispersion relation for 
graviton 

-fc%= (l + «2) (16) 
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3 Cosmological Dynamics and ^Ether Field 



3.1 Five-dimensional cosmological dynamics 

In this section we consider cosmological dynamics of 5-dimensional spacetime by applying 
Einstein general relativity to the product space, between 4-dimensional FRW-type spacetime 
and a circle S^. We assume the cosmological ansatz 

ds^ = -dt^ + aitfdx'dx^ij + hitfdy'^, (17) 

where i,j = 1,2,3, a{t) is the scale factor of non-compact 3-dimensional space, and h{t) 
denotes the radius of the compact fifth direction. The coordinates on are < ?/ < 27r. 
For our metric 0171) . the background solution for the equation of motion ([2]) can be written 
as 

m'^= (^0,0,0,0,^). (18) 

Using this background solution, the energy momentum tensor associated to the ^ther field 
defined in ([5]) can be written as 

■ iP' ■ b 1 

^"ol" = -y-^b' ^jl« = y-^fc'^V '^Wu = --^ + - ^HaHf, (19) 

We have defined the Hubble constants Ha = a/ a and Hj, = b/b, where dotted quantities 
represent the corresponding time derivative. As we mentioned in the previous section, Tab\u = 
when the extra dimension is stabilized 6 = 0. The fact that the ^ther field does not 
contribute to the energy density at the stabilized point implies that the ^ther field will not 
give any contribution to the effective potential of the radion. Hence other component such 
as Casimir energy is needed for stabilization of the extra dimension. However, as we shall 
see later on, the ^ther field can reduce the influence of the Casimir force. This property is 
important for stabilization mechanism when non-relativistic matter is present. 

Let us assume that the total energy-momentum tensor Tab\totai is decomposed into 

Tab] total I u ~l~ Tab \GC + Tab\p- (20) 

The contribution from non- minimally coupling to gravity Tab\GC is defined in equation f|T2|) . 
The component T^\p = diag{—p,pa,Pa,Pa,Pb) represents contribution from Casimir energy 
[3]. Casimir energy density p plays the role of 5-dimensional cosmological constant. Pa = —p 
and pb = — ^Q^2nb) ~ ~P ~ ^^bP are the pressure density in non-compact and compact 
direction respectively. By substituting Tab\totai into the Einstein field equation, we get the 
5-dimensional cosmological equations of motion 

3Hl + 3HaH, = 87rG{p+^v^H^), (21) 
3^-3HaH, = -87rG'|p + p5- (l-2a>2^}, (22) 

3^ + 9HaHb = 8TTG\^p + 2pt,-3pa-2{l-2agyAy (23) 
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where A = {\ + ^HaH^). 



3.2 Dynamics in the radion picture 

Since we are interested in our observed universe, it is useful to analyze the cosmological 
dynamics by considering 4-dimensional effective field theory. The equations of motion ( 12T|) - 
( l23l) can be obtained by varying the 5-dimensional Einstein-Hilbert action 

(fx^[j^R - -V,,V-^ + agU-u'R,, - V{h)). (24) 

V{h) denotes the potential term in 5-dimensional spacetime. Note that we omit the Lagrange 
multiplier term. For simplicity, we will set = in this section and this will not affect 
our main results. Let us start with KK-dimensional reduction of the above action from 5 
to 4-dimensional spacetime. Then, in order to make the resulting effective action in the 
canonical form, we apply Weyl rescaling g^^E = ^g^lu {p,^ = 1, • • • , 3) and define the new 
time variable dtE = V^dt, aEits) = V^a{t); ^ = 2TxbMl /m^i. Note that nipi is the Planck 
mass in 4-dimensional spacetime defined via the relation m^^ = {2'Khmin)Ml where 6mm 
denotes the stabilized radius of extra dimension. Thus f2 = 1 at 6 = hmin- The effective 
action takes the form 

/{ TDj^ 1 1 777-^ 2 >/i67r tj, 

^'xv^{^i?^--^7rV,vI>V.vl>--^^^e" ^-.^ V,V^-U{^)], (25) 

where U{^) = 27TbQ~'^V (b) is the 4-dimensional effective potential. Here we define the radion 
field ^ = ■^^V^ln{b/bmin) and = V^s = V^Ms. By using the background solution in 
(fTSjl . the above 4-dimensional action can be rewritten as 

/r m'^ 1 ^ 

d'xV^{Y^RE - -(1 + a')g^E'V,^V,^ - U{^)}, (26) 

where we define the dimensionless parameter = ^§f^- This action gives rise to the 
following set of equations: 

d^^ d'^ _ 1 dU 

d^+^^^d^ - ~jrT^)W ^^^^ 

Note that He = {daE/ dtE) / (ie is the Hubble constant in the Einstein frame. As we will 
explain later, the factor 1/(1 -f- a^) in the right-handed side of equation (|28l) weakens the 
effect of the potential gradient —dU/ and it is crucial for stabilization mechanism of the 
radion field ^ . To make contact with previous section, we note that the energy-momentum 
tensor associated with 5-dimensional action in fl2^ gives the relations 
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4 ^Ether Field and Casimir Energy 



We will start this section by reviewing the mathematical formulation to determine the 
Casimir energy for a scalar field, Ecas, and then investigating the effect of aether coupling 
to the Casimir energy. First, we consider Casimir energy of a non- interacting scalar field of 
mass, m, in M^^" x spacetime by following [151 [16]. We keep the number of non-compact 
spatial directions to be n for the moment and will set n = 3 at the end of our calculation. 
This scalar field obeys the free Klein-Gordon equation, 

(dad'' - m2)0 = 0. (30) 

The scalar field satisfies the periodic boundary condition in the compact direction, (f){y = 
0) = (f>{y = 27r). Its associated dispersion relation can be written as 

- k^k^ = m^ + —, (31) 

where, n G Z is the momentum number in the compact direction. Then, the total vacuum 
energy contributing to Casimir energy can be written as 



where = is the spatial volume of non-compact spacetime. Using the fact that 
/ f{k)d''k = 27r"/7r(n/2) / k''-^f{k)dk, we obtain 



where we define s = —{n + l)/2. Let us consider the massless case, m = 0. By using the 
zeta function regularization procedure, the Casimir energy density per one bosonic degree 
of freedom for massless scalar field can be written as 

^"^^^'^'^ = ^b = \l[J^^ ^''b''~'rr''~'a-2s + 1), (35) 

where ( denotes the zeta function and we take 27Tb to be the volume of compact dimension. 
For the massive case, we apply the Chowla-Selberg zeta function [16] in our regularization 
procedure and obtain the Casimir energy density per one degree of freedom for the massive 
scalar field: 

00 

pTa™ = -2(27r6)2^-i(m6)(i"2^)/2 ^n(2^-i)/2ir(i_2s)/2(27r6mn), (36) 

71=1 
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where Ki,{x) is the modified Bessel function. The fermionic degrees of freedom will contribute 
to the Casimir energy density with the same expression except for an extra minus sign. 

Let us consider the case that a scalar field couples to an aether field with a coupling con- 
stant a0. In the previous section, we showed that interaction with the aether field transforms 
the usual dispersion relation (|3T1) into its modified version ([8]). Accordingly, the Casimir en- 
ergy will be written as 



where we rescale k and m in such a way that A;^ = (1 -|- and im? = (1 + a|)m'^. 

By comparing Ecas{ci<j,) with the non-interacting Casimir energy Ecas, we see that the aether 
coupling rescales the Casimir energy and scalar mass by factors ^^^^ (i+c^2)-i/2 

respectively. Thus, we can immediately write down the Casimir energy density per one 
bosonic degrees of freedom with the aether coupling as 

rf-2s + 1)2^2^-1-3.-1 



7r 

r(-i/2) (1 + ^ 



massive i 
rboson \ 



for contributions from massless and massive scalar fields respectively. Other bosonic degrees 
of freedom contribute to the total Casimir energy in a similar way but with different coupling 
constants and masses. For example, the graviton in (n + 2)-dimensional spacetime has 
\{n + 2){n — 1) bosonic degrees of freedom with the ^ther coupling ag and mass m = 0. 
Using the modified dispersion relation for graviton f|T6l) . we can show that the graviton 
contributes ("'+2K"'~^) p^'^^'^^^ (ag) to the total Casimir energy density. Recall that, apart 
from modification of graviton's Casimir energy, non-minimal coupling of the aether field 
to gravity can also affect the dynamical evolution through the energy momentum tensor 

Tab I (GC) • 

For fermion case, we use the modified dispersion relation for the fermionic field in ffTOj) . 
We can show that the associated Casmir energy density per one fermionic degree of freedom, 
for both massless and massive case, is in the similar form of those from bosinic degrees of 
freedom with the over all minus sign and (l + a^) — > (1 -l-a^)^. The Casimir energy densities 
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for one degree of freedom of massless and massive fermion can be written respectively as 



Y[—2s + 1) 2^*6^'^~-'^7r^*^-^ 

P7ermion{a^) = T{-l/2) (1 + a^)2s C(-2g + 1), (40) 

The total Casimir energy density can be written in terms of sum over all degrees of freedom: 

AT ^massless/- \ , a/- rnassless f ^ \ , iyj rnassivef^ \ , aV rnassive f ^ \ I An\ 

P = J^bPboson + J^fP fermion ("vJ + ^bPboson l"^) + f P fermion{<^4^) ^ (42) 

where iVf, {Nf) and A*";, {Nf) are the numbers of bosonic (fermionic) degrees of freedom for 
massless and massive fields respectively. The nature of the total Casimir energy density 
depends on the relative magnitude of Nh, Nf, Ni, and Nf. 

In our model, N^ ^ 5 since, at least, the graviton is always present and it has five physical 
degrees of freedom in five-dimensional spacetime (n = 3). The compact fifth direction would 
not be stable if there is only the graviton field in the bulk, i.e. it will collapse to Planck size, 
due to the negative Casimir energy associated with quantum fluctuations of the gravitation 
fields. Therefore, it is natural to add more positive contribution to the Casimir energy by 
assuming that there are fermions in the bulk. However, we cannot create the minimum 
of p by including only the massless fermionic fields i.e. the Casimir force is attractive for 
Nf < Nb and becomes repulsive when Nf > Nb. Hence Nf = Nb does not give us any stable 
fixed point. We should expect a minimum to be produced if we include massive fermionic 
degree of freedom [T7] . This can be explained qualitatively as the following. Let us consider 
Casimir energy of a fermion with mass M. For the region where b -C 1/M, the vacuum 
energy should have the same form as in the massless case. In particular, for Nf > Nb, the 
net Casimir force will be repulsive. In the other region where b ^ 1/M, the contribution 
from the massive fermion mode is negligible compared to that of the graviton and the total 
Casimir force becomes attractive. Hence, there must be a stable fixed point between these 
two regions. 

In this paper, we consider a toy model where the particle spectrum in the bulk consists 
of a bulk graviton, a massless Dirac fermion, a massive Dirac fermion with mass rrif, and 
eight massive scalars with equal masses rris = Xrrif. Here A is the mass ratio. The presence 
of the massless fermion and the massive scalars is to ensure that the vacuum energy at the 
minimum has positive value, i.e. pmin > 0. We summarize the particle content in the bulk 
in Table 1. 

There is no unique choice of bulk particle spectrum for this purpose, the other combina- 
tions of bulk fields, for example in Ref. [3], can probably create minimum for the vacuum 
energy. Our particular choice is convenient for investigating the effect of the aether-matter 
interactions on the Casimir energy of the bulk fields. Note that we do not attempt to justify 
the existence of these bulk fields phenomenologically because we want to demonstrate that 
stability of the extra dimension could be achieved, if these particles are present in the bulk. 
The Dirac fermions in five dimensional spacetime have eight physical degrees of freedom. For 
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Particles 


Degrees of freedom 


Mass 


Coupling constant 


a bulk graviton field 
a massless bulk fermion field 
a massive bulk fermion field 
8 massive bulk scalar fields 


Nb = 5 
Nf = 8 
Nf = 8 
Nb = 8 





nif 
iTis = XrUf 


a, = (16vr^)V2 



Table 1: The particle spectrum in the bulk, their degrees of freedom, their mass and the 
coupling constants characterizing their interaction with the ^ther field. For simplicity, we 
assume the universal fermionic coupling for both massless and massive Dirac fermion. 



more realistic models which have chiral fermion on the brane, we can impose the orbifold 
reflection symmetry: y —y on the compact direction. However, this will not affect our 
main results on stabilization of the extra dimension. We will ignore this issue for simplicity. 

5 Effects of The ^Ether Field on Stabilization of The 
Extra Dimension 

5.1 Stabilization in vacuum dominated universe 

We first consider the universe where there is no non-relativistic matter and the Casimir 
energy density is the dominant contribution. Let us start with the case where there is no 
interaction between aether field and bulk matter i.e. setting a^i, = = ag = 0. The plot 
of the total Casimir energy density in 5-dimensional spacetime p using the expression in 
equation fH2|) and particle spectrum in Table 1, with the mass ratio A = 0.516, and its 
corresponding 4-dimensional effective Ucas = 27r(6^j„/6(\E'))p(^') is given in Figure [H The 
local minimum of Casimir energy density, pmin = 23.4316(m//40)^ is located at bmin = 
0.01461(40/mj). In order to get the positive minimum and the radion can stabilize, we 
must choose the value of A from a very narrow range, 0.516 < A < 0.527. By solving the 
5-dimensional equations of motion (!2Tl) -( |23l) numerically, we can demonstrate that the extra 
dimension is stabilized at the radius bmin as shown in Figure [2l Notice that we set the 

expansion time scale to be in the unit of Hubble time = H^q = ^?>rn^j8Txpc ~ 10^° 
years. The critical density pc can be written in terms of the minimum value of Casimir 
energy density pmin as pc = (1 + 0.24/0.76)(27r6^i„)p^i„. 

By comparing 4-dimensional effective Casimir energy density p[^|„ = i^nhmin) Pmin with 
the observed value of energy density for dark energy, p^^|, ^ (2.3 x 10^^ eV^)^, we get mf ^ 
4.18 X 10~^ eV. Then, the radius of extra dimension bmin ~ 13.96 eV~^ ~ 2.75 x 10~^ m. 
This leads to the quantum gravity scale in the bulk, M^: ^ 1.19 x 10^ GeV. Note that we do 
not attempt to address the mass hierarchy problem in this paper. In order to compare with 
the ADD brane world scenario, it is better to be demonstrated with 6-dimensional models 
as shown in [3j. Since our aim is to study the role of aether fields on stabilization of the extra 
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Figure 1: Left: Casimir enery density p is presented in the y-axis in units of (m//40)^. The 
X-axis is the radius of the fifth dimension h in units of 40/mj. Right: The 4- dimensional 
effective Casimir potential Ucas = 2vr(6^j„/6(\E'))p(\E') in unit of (m//40)*^ is presented in the 
y-axis. The x-axis denotes in the unit of nipi. Here we set a = = = = 0. 

dimension, the 5-dimensional model is good enough for our purpose. We will leave the mass 
hierarchy problem for future works. 

The role of aether field on dynamical evolution of the extra dimension is also illustrated 
in Figure [21 We can see that, as the value of v increase, the moduli field b feel less potential 
gradient. Its oscillation frequency and amplitude decrease. In our caculation v is in the unit 
of (mj/40)^/^t//. From the previous paragraph mj = 4.18 x 10~^ eV, this gives 1 unit of v 
is equivalent to (1.02 x 10^ GeVY^"^ ^ (0.09M=^)^/^. At very high v, for example v = 100 
or approximately ~ (9M*)^/^, the period of oscillation is so long that b reaches equilibrium 
before showing any oscillating behavior. The scale factor b rolls slowly to its stable fixed 
point. 

Let us consider the situation where aether field couples to the bulk matters. As we 
discussed earlier, interactions with the aether field reduce the effective mass of the bulk 
fields. For example, the effective mass for scalar field of mass rris would be 

^'(e//) =^s(l + «J)"'- (43) 

This will alter the shape of the potential as we demonstrate in Figure [31 Interestingly, there 
is an advantage of coupling the bulk fields with the aether field. It seems that we get the 
wider range of parameter space for the mass ratio A that allows positive minima Pmm > 0, 
i.e 0.05 < A < 0.80. 

5.2 Stabilization in the universe with non-relativistic matter 

In this section we consider the role of aether field on stabilizing mechanism of the extra 
dimension in the more realistic model of our universe i.e. a model containing non-relativistic 
matter. Let us fist demonstrate the destabilizing effect due to non-relativistic matter by 
following Greene and Levin in Ref [3]. We assume that there is matter living in the bulk. 
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Figure 2: The dynamics of the scale factor b{t) for the compact direction as a function of 
time with different values of parameter v. In the absence of aether field v = 0, b{t) shows 
oscillation behavior around its critical value bmin before stabilizing at this value. Non-vanish 
value of V reduces the influence of Casimir force. As the value of v increases, the oscillation 
frequency and amplitude decrease. If the vev of the aether field is large enough, for example 
V = 100, oscillation behavior disappears. The extra dimension evolves smoothly to its stable 
fixed point. The time variable t is presented in the unit of Hubble time tn- The time for 
stabilization to occur is around ~ Gtn- 



The condition for stabilization of 6 is 17 ^ 10 ^. 
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Figure 3: Interactions between bosons/fermions and aether field can affect the Casimir energy. 
In this figure, we fix the mass ratio A = 0.516. The Casimir enery density p is presented 
in the y-axis in units of (mj/40)^. The solid, long dashed and short dashed line denote the 
Casimir energy density when the coupling constant (a^,^^) = (0.0, 0.0), (1.0, 0.644), and 
(1.5, 0.897) respectively. The value of bmin and pmin increase as we increase the value of the 
coupling constants. The shape of the potential well gets shallower as the coupling increases. 
We set dig = for simplicity. 

This can be done by adding a 5-dimensional matter term into the 5-dimensional action in 
equation (1211) 



This is equivalent to adding the matter energy density pm oc l/2Tio?h into the 5-dimensional 
cosmological equations of motion The matter energy density pm includes contri- 

bution from baryonic matter and cold dark matter. By comparing with the observational 
data and supposing that all dark matter is cold, the matter density today p^o is roughly 
26% of the total energy density of the universe. The Casimir energy density will be respon- 
sible for the other 74% of the total energy density today in the form of dark energy, (pao)- 
Thus, we have the relation, p^o = (2.6/7.4)pAo- Note that the energy density of dark energy 
in our 4-dimensional observed universe today can be written in terms of the minimum of 
5-dimensional Casimir energy density and the stabilized radius of the extra dimension as 
Pao = Pmm(27r6mm) = (2.3 X lQ~^eVY. By using (ao/a) = 1 + 2;, oq is the scale factor today 
and z is the red-shift, we get 




(44) 




(45) 
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In this case, the 5-dimensional equations of motion fl21l) - fl23l) become 



ml + maH, = 87rG(p + p„ + it;X), (46) 
3^-3HaHb = -8nG\^p + p^+pb-il-2agyAy (47) 

3^ + 9HaH, = 87rG|p + p^ + 2pb-3p,-2(l-2a>2^}. (48) 

From equation (HHl) . the stabihty conditions {Hb = 0, Hb = and A = 0) require 

Pb = -2p - -pm. (49) 

By using equation fHTl) and requiring that f > when h is stabihzed, we get the constraint 
Pm < 2p. This is the same constraint that we have for a (3 + l)-dimensional vacuum 
dominated universe. Thus, as pointed out in [3], this model describes the (3 + l)-dimensional 
vacuum dominated universe when the extra dimension is stabihzed. 

By using the conservasion of the energy-momentum tensor in four and five dimensions, we 
can easily show that the radion field will be driven toward the minimum of the 4-dimensional 
effective potential 

2 (4) /, \ 2 

TT TT I Pl^ Pm jr , Pm Omin \ 

Ueff = Ucas + Y = ^Cas + -r\^] > (50) 



n 4 ' 4 . 

where we define the 4-dimensional matter density pm = Pmi^Tfb) = j^Pmin {'^'^i'min) (1 + 
which is a function of (1 + z)^ and does not depend on the radius of the extra dimension h. 

Numerical results for the effective potential Ueff are illustrated in Figure HI Here we 
choose A = 0.516, and ignore the interaction terms by setting = = ag = 0. At 
^ = 0, the presence of non-relativistic matter would lift up the minimum of Ueff slightly. 
However, at early time, high red-shift, the 1/6^-term in (1501) becomes dominant and destroys 
the presence of the minimum. This effect will drive h to expand even though there is a 
minimum today since the radion field ^ has already rolled pass the minimum and cannot 
get back to the stable point. Notice that this effect is the same if matter is confined to the 
brane. This can be shown by adding a 4-dimensional term into the 5-dimensional action 

S = Sr,i:) + j d'^X^/^g^Cmatter- (51) 

There is no dimension reduction in the matter term but the conformal transformation rescale 
the matter energy density as p^ /Vt^. Thus, the total effective potential is the same as in 

dSoD. 

There is an interesting phenomena if we turn on the ^ther field. By carefully tuning 
its norm, we can show that h can be stabilized even if there is non-relativistic matter in 
the universe. This can be shown by solving equations of motion ( I46l) -( l48l) numerically. The 
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Figure 4: This is the plot of the effective potential f/e//(^) at red shift z = 0.0, 7.0, 9.0, 
11.0, 13.0 in the unit of (mj/40)^ and A = 0.516. The local minimum of Ueff{'$) no longer 
exists when the red shift is increased. 

result is shown in Figure O where we choose A = 0.516, rrif = 40 and f = 10 (or equivalent to 
V ~ (0.9M^,)^/^). We set the initial conditions such that b{t = 0) = 0.606mm and ignore the 
effects of interaction terms by setting a^j, = = Og = 0. From Figure HI the minimum of 
Ueff{^) disappear when z > 11. So, in order to make sure that the minimum of Ueff{'^) does 
not exist in the initial configuration, we choose the initial time to be at z = 10, 000 which is 
approximately the time of matter-radiation equality (teg)- Moreover since the universe also 
contains radiation component which is not considered in this model, it is reasonable to set 
tinitiai = teq = 0. However, if we ignore the radiation part, we can go back to the much earher 
time and show that there are some values of v that the extra dimension can be stabilized. 

The role of the aether field can be explained qualitatively in the radion picture. The 
aether factor 1/(1 + a^) on the right-handed side of equation (!28|) reduces the influence of the 
potential gradient —dUeff/d'^. As a consequence, it will slow down the oscillation frequency 
of \E' around the minimum of the potential Ueff{^)- If this factor is big enough, ^ will move 
down the potential at very slow speed. We can tune v such that there is enough time for 
the universe to create the minimum of Ueff{^) before the radion rolls pass it. By this 
mechanism, stability of the extra dimension can be restored. 

Let us compare the stabilization time tgtab of the moduli field with the age of the universe. 
The age of the universe in our model is 



where we set ^casimir = 0.76 and VLm = 0.24. From Figure O the stabilization time tgtab ~ 
Qtu- Then, tgtab ~ 3.90tage is greater than the age of the universe. 




(52) 
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Figure 5: These graphs illustrate cosmo logical dynamics of the universe which includes non- 
relativistic matter content and ^ther field in the extra dimension. Left: The scale factor a 
(upper) and the Hubble constant for the non-compact dimensions Ha (lower) as the function 
of time. The fact that Ha oscillates indicating the period of deceleration and acceleration 
before it settles down to the constant value and a enters a de Sitter phase. Right: The scale 
factor h (upper) and the Hubble constant for the compact extra dimension H^ (lower) as the 
function of time. H^ oscillates between positive and negative region before settles down to 
zero. The extra dimension is stabilized although non-relativistic matter is present. 
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6 Conclusions and Discussions 



In our 5-dimensional model, we have shown that the aether field reduces the infiuence of the 
potential gradient and slows down the oscillation frequency of the compact extra dimension. 
For vacuum dominated universe, the Casimir energy from the extra dimension acts as a 
stabilizing potential for the moduli field while driving accelerated expansion in the non- 
compact directions. The ^ther field will slow down the oscillation behavior of the moduli or 
even smooth it out. For the universe which non-relativistic matter is present, this effect is 
proved to be crucial for stabilization of the extra dimension. If the vev of aether field is of 
the order of the 5-dimensional Planck mass, v ~ 0{mI'^), it can slow down the evolution 
of the moduli field such that there is enough time to create the minimum for the effective 
potential. 

In this paper we assume homogeneous and isotropic distribution of non-relativistic mat- 
ter. However, local matter distributions might perturb the radion and knock it over the 
minimum, causing the (local) catastrophic expansion of the fifth dimension. In [3], it was 
also noted that the minimum of the potential well is generally not deep enough to prevent 
the quantum tunneling of the radion. At this stage, it is not clear whether these two diffi- 
culties can be solved by the new mechanism. These aspects of instability in the presence of 
the aether are still open questions. 

Note that the constancy of the 4-dimensional cosmological constant up to very early 
epoch of the universe will post strong constraint on the size of the extra dimension. The 
oscillation behavior of the moduli field may contradict with astronomical observations. In 
order to construct a more realistic cosmological model of this scenario, the extra dimension 
should reach its stable fixed point before the present time i.e. tstab ^ iage which require fine 
tuning of many parameters. The new possible solution is that we assume very high value 
of V so that the oscillation of h has a very long period. The moduli will evolve smoothly 
with no oscillating behavior. We can choose the value of v such that the size of the extra 
dimension changes so slowly and it cannot alter the results of the Big Bang model. 

Another interesting idea is to imagine that the universe started with a very symmetric 
state which all spatial dimensions are compactified with the equal radius, for example, topo- 
logically an 4-dimensional torus. In general, the Casimir energies in this compact universe 
generate stabilized potential for the radius of all directions. On the other hand, at the very 
early time, the energy density of matter and radiation will be the dominant contribution. 
This will destabilize the moduli fields and all directions will become large. However, if the 
Lorentz symmetry is spontaneously broken in some direction i.e. there is a non-vanishing 
aether field pointing in the fifth direction, it will slow down the dynamics of moduli field 
associated to the broken direction. The broken direction will be compactified at stabilized 
radius while the unbroken directions are allowed to expand. This cosmological scenario may 
establish a connection between the dimensionality of spacetime and the violation of Lorentz 
symmetry. We leave this issue for future investigation. 
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